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Invariant tori play a fundamental role in the dynamics of symplectic and volume-preserving 
maps. Codimension-one tori are particularly important as they form barriers to transport. Such 
tori foliate the phase space of integrable, volume-preserving maps with one action and d angles, 
i—i For the area-preserving case, Greene's residue criterion is often used to predict the destruction 

1 ^ i of tori from the properties of nearby periodic orbits. Even though KAM theory applies to 

the three-dimensional case, the robustness of tori in such systems is still poorly understood. 
We study a three-dimensional, reversible, volume-preserving analogue of Chirikov's standard 
map with one action and two angles. We investigate the preservation and destruction of tori 
under perturbation by computing the "residue" of nearby periodic orbits. We find tori with 
Diophantine rotation vectors in the "spiral mean" cubic algebraic field. The residue is used to 
\Q generate the critical function of the map and find a candidate for the most robust torus. 

in 

^ 1 Introduction 

This paper is concerned with the breakup of invariant tori for volume-preserving maps. Tori are 
prominent as invariant sets for many dynamical systems. In particular a dynamical system is 
deemed to be "integrable" when almost all of its orbits lie on tori and the dynamics on each 
torus is conjugate to a rigid rotation with some rotation vector uj. This is the Arnold-Liouville 
formulation of integrability |Arn78] . and as such holds for Hamiltonian flows and symplectic maps, 
but also is essentially the notion of "broad integrability" introduced by Bogoyavlenskij for more 
general systems |Bog98|. 
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Here we will consider families of maps f £ on M = T d x R k of the form 

x = x + Q(z) — eh(x, z, e), 



(1) 

z = z — eg{x, z, e), 



that have an "integrable limit" when e = 0. Taking T d = R"/Z , the x coordinates represent 
d-angles with unit period, and the z-coordinates represent k action-like variables. Thus the map 
([!]) is said to be a d-angle, /c-action map. We will assume that / is volume-preserving, i.e., that its 
Jacobian satisfies 

detDf £ = l. 



Examples of such maps include Chirikov's standard map (with d = k = 1) 

x' = x + z — esin(27ra;), 
z = z — esin(27rx). 



(2) 



This map is area-preserving, and is thus both volume-preserving and symplectic. Another oft- 
studied case is the four-dimensional symplectic map due to Froeshle (with d = k = 2) [FS73]. 
The dynamics for d + k = 3 applies to the motion of a passive scalar in an incompressible fluid 
[PF881 ICFP96] , as well as granular mixing |MLO07| and magnetic field-line flows [TH85J . In this 
paper we will study a two-angle, one-action map — a (2 + l)-dimensional map — that models the 
typical dynamics near a rank-one resonance |DM12j . Our model is introduced in ^3j 
The frequency map 

Q : R k -»• R d (3) 

plays a key role in the dynamics of f £ . In particular every orbit of the integrable limit, /o, lies 
on a d-dimensional torus, T z = T d x {z}, on which the dynamics is simply a rigid translation, 
x t = %o + with uj = £1(zq). We call any ti-torus that is homotopic to 7o, rotational. When the 
rotation vector uj is incommensurate, namely 

{( P ,q)€Z d xZ:p-u; = q} = {(0,0)}, 

such orbits are dense on T z . KAM theory studies the persistence of the rotational tori of perturba- 
tions of integrable system like /o [BHTB90, Pos0l| lDLL01| . The application of this theory requires 
that f e preserve some structure (e.g., is symplectic or reversible), is smooth enough, has a nonde- 
generate frequency map, and finally that the rotation vector of the tori be sufficiently irrational, 
e.g., to G D s is Diophantine. 

V s = Lo G M. d : 3c> s.t. \p ■ ui - q\ > ^\ , (4) 

for some s > d. Since we assume that k = 1, whenever a rotational invariant torus exists it has 
codimension-one and therefore is a barrier to transport. 

The orbits of fo with rational rotation vectors, to € Q d , are periodic. More generally, each 
periodic orbit of ([!]) can be assigned a rotation vector by lifting the map to R d x R fc . The lift of a 
period-n orbit obeys 

(x n ,Zn) = (x + m,z ), (5) 



2 



for m £ Z d , and thus has rotation vector uj = m/n. We call these (m, n)-periodic orbits, see £3.1 

Following John Greene |Gre79| , we study the persistence of rotational tori by finding a sequence 
of rational vectors mg/ne that limit on a given incommensurate vector uj. Greene hypothesized that 
the limit of the corresponding (me, n£)-periodic orbits of f £ is a quasiperiodic orbit dense on an 
invariant torus if and only if the periodic orbits remain stable in the limit; more specifically, if their 
residues remain bounded. This is Greene's residue criterion, see ^2] As we recall in that section, 
proofs of this result are known in certain situations. 

Moreover, Tompaidis generalized Greene's criterion to higher-dimensional symplectic and quasi- 
periodically-forced symplectic maps [Tom96aj. Indeed, the residue criterion has been applied to 
study the breakup of two-tori for several four-dimensional models, e.g., the Froeshle map [BM93, 
IBS931 ITom96bl \KA97\ ITom991 IZHSOlal IZHSOlbl IGFL04 ] , and a quadratic map |VIB97j . We are 
not aware of previous use of the residue criterion for volume-preserving maps apart from a quasi- 
periodically-forced, area-preserving map [AC S91| IACS92} ITom96b] . a (2 + l)-dimensional map with 
one of the components of f2 set to a fixed irrational value. 

As we discuss in £3.2, our model — like the standard map ^ — is reversible, and this makes 
finding periodic orbits especially easy: a two-dimensional secant method suffices, see §|4j In addition, 
since one of the multipliers is always 1 for symmetric orbits, there is a natural definition of the 
residue, see §3.4| 

One complicating feature of ([I]) with k < d is that the image of the frequency map Q is at most 
a fc-dimensional subset of the <i-dimensional space of rotation vectors; a consequence is that orbits 
with given rotation vectors typically do not exist for fixed parameters. In our (2 + l)-dimensional 
model, we get around this by adding a new parameter, 5, to Q so that : (z,5) i-> (ui\,uj2) is a 
diffeomorphism. With this modification, we observe, in Q that there are symmetric periodic orbits 
for any e and any m/n G Q 2 . 

We will show in ^5] that the residues of high period orbits appear to undergo a rapid transition 
from "nearly zero" to "exponentially large" as e grows, just like for the standard map. In ^6] we 
use this transition to get a reasonably sharp estimate for a critical set of parameters at which a 
given torus is apparently destroyed. 

The tori that we study in ^6] have rotation vectors that are integral bases for a cubic algebraic 
field Q(cr), where a is the "spiral mean" [K086j . see £6.3 It has long been conjectured that simple 
cubic irrationals like the spiral mean could be the analogue of the golden mean for the standard 
map Q. Indeed, most of the previous studies of Greene's criteria for multidimensional tori have 
used vectors in Q(a), though there have been several that studied other cubics |Tom96bt IZHSOla) 
IZHS01b| and even quartic irrationals [V1B97L lKA"97j . 

Another remarkable conjecture in |Gre79j is that the last invariant circle of ^ has the golden 
mean rotation number. There is strong numerical support for this conjecture, and more generally for 
the conjecture that circles with "noble" rotation vectors appear to be locally most robust [MS92], 
Do these notions have a higher-dimensional generalization? We are not aware of any previous 
progress on this question. In £j8]we look for the last torus for our model by finding the set of critical 
parameter values, e CT (u)), for a set of tori whose rotation numbers are integral bases of Q(cr). We 
use Kim and Ostlund's generalization of the Farey tree, see £6.1, to systematically generate sets of 
Diophantine rotation vectors. 
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2 Greene's Residue Criterion 



John Greene studied the persistence and destruction of rotational invariant circles of the standard 
map Q |Gre681 IGre79| by approximating them with sequences of periodic orbits. The Poincare- 
Birkhoff theorem implies that the standard map has at least two (m, n)-periodic orbits for any 
choice of rational rotation number — |Mei92j . These orbits, which exist for all e, are often called 
the Birkhoff orbits. Since the standard map is area-preserving, the product of the multipliers of 
any period-n orbit is 1, and so its stability can be completely characterized by a quantity Greene 
called the residue: 

R = |(2 -r), where r = tr(D f n (x, z)). (6) 

The residue conveniently encodes the stability of an orbit: it is elliptic (complex, unit modulus 
multipliers) when < R < 1, hyperbolic (real, positive multipliers) when R < 0, and reflection 
hyperbolic (real, negative multipliers) when R > 1. One of the two Birkhoff orbits of ^ has 
positive residue while the other has negative residue, and Greene showed that for a period-n orbit 
of (§, 

R = (D(e n ), (7) 

both for £ <C 1 and for e> 1. 

Greene's 6 th assertion, now known as Greene's Residue Criterion, is perhaps the most aston- 
ishing of the conjectures in |Gre79j . In particular, consider an irrational rotation number oj G M, 
with continued fraction u) = [ko, k\, ■ . .] = fco + l/(&i + . . .), an d let 

— = [fco,*i,...,fci] (8) 

ni 

be its I th convergent. Thus the (mg, n£)-Birkhoff orbits, with residues R» , have rotation numbers 
that converge to u. One formulation of Greene's criterion is: 

Greene's Residue Criterion: A rotational invariant circle of an area-preserving twist 
map with a given irrational rotation number oj exists if and only if the residues of its 
convergent Birkhoff orbits, , remain bounded as — — > oj. 

A stronger version of this criterion asserts that when there is an invariant circle the mean residue, 

fi((j) = lim -log |^|, (9) 

l— >oo ni 

exists and is negative for any sequence of (mg, n£)-orbits whose rotation number converges to oj. 
Greene studied in particular the invariant circle with golden mean rotation number 

(t> = l(i + Vh) = [1,1,1,...]. (io) 

For the convergents to <fi, numerical studies show that there is a parameter value £ C r(<fi) such 
that Rf -)■ as I — > oo whenever e < e cr (4>). Conversely, whenever e > e cr ((p) the residues grow 
exponentially with the period, and there appears to be no golden invariant circle. Greene estimated 
e cr by computing a sequence of parameter values for which \R^\ reaches some fixed value Rth > 0, 
obtaining 

e cr (<P) pa - 1 -0.97163540631. (11) 

Z7T 
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The value of the threshold Rth is irrelevant, but Greene found that Rth ~ 0.25 gave the most rapid 
convergence. 

Some aspects of Greene's residue conjecture have been proven. 

Theorem 1 (Residue Criterion for Twist Maps [Mac92, FdlL92j). Suppose f £ is an analytic, area- 
preserving twist map, and the sequence of (me, ni) -orbits converges to an analytic invariant circle 
on which the dynamics is analytically conjugate to rigid rotation with rotation number u G T> s , then 
there are constants C, K > such that 

\Ri\ < Cexp (-K\u - mf/^l- 1 /' 1 ^)) 

In particular since \oj — m£/n£\ < c/nj for continued fraction convergents, this implies that Ri — > 
exponentially in the period. The residue criterion also applies to twist-reversing maps [DdlLOO], as 
studied numerically by |dCNGM96l lAWM03j . 

Aspects of the converse statement of the residue criterion have also been proven for area- 
preserving twist maps. From Aubry-Mather theory, these maps have a set of "minimizing" orbits 
for each uj, and when oj is irrational this set is either a circle or a Cantor set — a cantorus. If the 
cantorus has a positive Lyapunov exponent, then there exists a sequence {m^ng) such that the 
mean residue ^ is positive and has the value of this exponent [FdlL92j . Moreover, if the cantorus 
is uniformly hyperbolic, then this sequence can be taken to be the sequence of minimizing orbits 
(the negative residue Birkhoff orbits) |Mac92| . 

The residue criterion has also been generalized to higher dimensions. For example, a 2d dimen- 
sional symplectic map has partial residues 

R(j) = l(2-X j -Xj l ), j = l,...,d 

for each of the reciprocal pairs of multipliers of a given orbit [KM89]. Tompaidis [Tom96l] proved 
that if f £ is a C T symplectic map with r > 1, and the twist, DO, is nondegenerate then, when 
there is a C r , Diophantine, invariant d-torus, the partial residues of any (m, n)-periodic orbits with 
rotation numbers sufficiently close to oj obey the bound 

\R^\ < Cn\nu - m\ k 

for any positive integer k < (r — l)/2s. His results also apply to the rotating standard map of 
|ACS92j . 

We are not aware of any theorem, however, that applies to volume-preserving maps more 
generally. Thus we turn to numerical investigations of a model. 



3 One-Action Maps: A Three-Dimensional Model 

We will study the rotational tori of a (2 + l)-dimensional map of the form 0. When k = 1 the 
image of a frequency map Q for is a curve. Near resonance, p-Q(z) = q, the rotational invariant 
tori of ([T]) are typically fragile. The resulting local dynamics is strongly influenced by whether the 
curve tt(z) is transverse to or tangent to such a resonance [DM12]. We will study dynamics of the 
two-angle case with the frequency map 

n(z,5) = (z + 1 ,f3z 2 -5), (12) 
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which, as was shown in [DM12J, is a normal form for dynamics near a "rank-one" resonance. When 



j3 7^ 0, the map (12) satisfies the nondegeneracy condition 



det(D z n,D 2 z n,...,Din) 



(13) 



that is sufficient for KAM theory [CS90a, Xia92]. For fixed (7, /3, <5), the image O : z 1— > ui is a 
parabola. However, when thought of as a map f2 : (z,5) 1— > (^1,^2)1 ^ becomes bijective. The 
parameter 5 is also important in another sense: it unfolds a tangency with the (0, 1,0) resonance 
(at z = 5 = 0). As was shown in [D M12] , near a resonant tangency, even though the "twist" 



condition (13) applies, the dynamics is like that of nontwist maps. Consequently, we will take 5 to 



be an essential parameter, but will fix j3 and 7. 

Another requirement for the preservation of the rotational tori of ([I]) for nonzero e is that the 
map satisfy an "intersection property" . A necessary condition is that the force has zero average, 



g(x, z, e)dx = 0. 



(14) 



Indeed, if (14) is not satisfied, then f e may have no invariant tori; for example, when g(x,z,e) 



const 7^ 0, then the z coordinates of ([!]) drift and there are no recurrent orbits for any e / 0. For 
the two-angle case, we will use the simple form for the force, 



g{x) = asin(27rxi) + 6sin(27TX2) + csin(27r(xi — £2)). 



(15) 



The three terms in g represent resonant forcing for (p, q) = (1, 0, q), (0, 1, q) and (1, —1, q), respec- 
tively, for each g£Z. The last term explicitly couples the two angles (xi,X2); however, note that 
even when c = these are coupled through the frequency map ( fl~2] ). 

As a final simplification, we will choose the perturbation h in ([!]) so that eh = Vt{z) — fl(z'), 
giving the model of [D M12] 

x = x + Q(z', 5) modi, , , 

z = z — eg{x). 

An advantage of this form is that it is always a homeomorphism; indeed, it has the inverse 

f-\x,z) = (x-n(z,S),z + eg(x-n(z,S))). (17) 

In addition, this map is an exact-volume-preserving diffeomorphism whenever Q and g are C 1 and 
g satisfies (14). 

One could also regard / as a diffeomorphism on T 2 x M 2 , by adjoining the trivial dynamics 



5' = 5 to (16). 



Following |Meil2j . we will typically use the "standard" set of parameters 

a = b = c = l, (3 = 2, 7 = i(\/5-l), 



(18) 



and think of e and 5 as the parameters that govern the strength of the forcing and the unfolding 
of the tangency, respectively. 



Unless otherwise mentioned, all the computations below are for (16) with the frequency map 



(12), the force (15), and the standard parameters (18). 
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3.1 Periodic Orbits 



Rotational periodic orbits of (16) are partially classified by their rotation vectors. Let F be the lift 
of / to the universal cover M 3 of M obtained by simply removing the mod 1 from (16). A sequence 



r = {{x t , z t ) = F(x t -i,z t -i) :t£Z} (19) 

is a type (m, re) £ Z 2 xN periodic orbit of F if it obeys Alternatively, noting that the rigid 
translation operator 

T m (x, z) = (x + m, z), (20) 

is a symmetry of F, then the periodicity condition becomes T_ m F n (x, z) = (x, z). If m and n are 
coprime, i.e., gcd{m\,m2, n) = 1, the projection of each such orbit onto M is a period-n orbit of /. 
The rotation vector of an orbit is the average increase in the angle per iteration, 

w(r) = lim -(xt - xq) 

t— >oa t 



if this limit exists. Each (m, n)-orbit has rotation vector 



An (m, n)-orbit of (|16[) with the frequency map (12) must satisfy the three equations 



n-1 



n = y^n(zt,6) 

t=o 
n-1 



mi 
m 2 



i 1 

n i + E t =o 

-n5 + pJ2 



zt 

n-1 9 
t=0 z t 



(21) 



t=o 



with (xt,zt) = F t (xo,zo)- Ideally a system of three equations in three unknowns has isolated 
solutions; however, we recall again that for any fixed 5 and a given (m,n), there are typically no 
solutions of (21). For example, when e = 0, the curve Q(z,S) will typically not intersect the point 



in 
n 

for e 



and there will be no (m, n)-orbit. However, including the parameter 5 in (12), there is an orbit 
at 

[z ,6 )- 



m 2 \ 



( mi a 

l,pz 

V n n / 

and any value for x. This corresponds to a two-torus of (m, n)-orbits of / 



(22) 



More generally, if an orbit satisfies the first and last equations in (21), then the value of 5 is 



completely determined by the second of these equations, which can be rewritten as 



n 



n-1 

t=o 



zf — 0J2 = (z 2 ) — 0J2- 



(23) 



Thus 6 is fixed by the mean square average of the action along the orbit. Of course the system (21 ) 



cannot be solved independently of 5, since the iterates of (16) depend on its value. As we will see 
below, with 5 included as a varying parameter, it appears that f e has orbits for all (m, n) £ Z 2 x N. 

For example, fixed points of ( 16) occur with (z* ,5*) given by (22) upon setting n = 1. Thus these 
values are uniquely determined for any (7711,7712). The angles are then determined by g{x*) = 0, 
which generically has a one-dimensional set of solutions that project to circles on T 2 . For example, 



7 



-0.4 0.0 0.4 



Figure 1: Curves of fixed points for (16 1 projected onto the angles x. The thick (blue) lines correspond to the 

(1.0,0.7,0.3). The stability of these orbits is indicated for 



case a — b = c, and the thinner (red) curves to (o, b, c 
e = 0.7 and mi = 1; unstable orbits with R < are shown as dashed lines and those with R > 1 as dotted lines 



if a = b = c fixed points occur on the three circles x\ = 0, x 2 = \, and x\ — x% - ,, 
there is at least a pair of circles of fixed points; an example is shown in Fig. [TJ 
For period two, the system (21) becomes 



\ . More generally 



Xl 

S 



leaving one transcendental equation 



|(mi + sg(x )) -7, 
|(mi - sg(x )) - 7, 

X + fi(2l, 6), 

\{l3zl + I3zl-m 2 ) 



ffOo) + g(xi) = o 



(24) 



for xq. Consequently, there is again a set of curves of period-two orbits. Note however, that the 
the value of 5 depends upon e and m for n = 2. This holds more generally for higher period orbits 
as well. 



3.2 Reversibility 



As for the standard map Q, it is very convenient that the map (16) with (15) is reversible as this 
simplifies the search for periodic orbits. Recall that a map is reversible if it is conjugate to its 
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n 


Si 




Sf 




X 


even 


Si 


Si 


oT_ 


m 


n 
2 




S 2 


s 2 


oT_ 


-m 


11 

2 


odd 


Si 


S2 


oT_ 


m 


n+1 
2 




S2 


Si 


oT_ 


m 


n-1 
2 



Table 1: Initial, Si, and final, Sf, symmetries for (m, n)-symmetric periodic orbits for a lift F of a reversible map 
with reversors Si and S2 = F o Si and discrete rotation symmetry T m . The initial point (zo, zo) £ Fix(Si) maps to 
the point (x x ,z x ) £ Fix(5/) in x iterations. 



inverse, that is if there exists a homeomorphism S such that f o S = S o f^ 1 . The collection of the 
symmetries and reversors of a map form a group, the reversing symmetry group |RQ92 |LR98] . 
Whenever g is an odd function (as in (15)), the map (16) is reversible, with the two reversors 



Si(x,z) 
S 2 (x,z) 



(-x,z- eg(x)), 
(—x + tl(z, 5), z) 



(25) 



To see this note that 



since the inverse of / is (17) 
consequently 



S 2 (x,z) = foS\(x,z) = S x of 1 {x,z), 

Since g is odd, both Si and S 2 are involutions, that is Sf 

f = S 2 o Si. 



id. 



Factorization into a pair of involutions also holds for Chirikov's standard map ([2]), and this was 
exploited by Greene and MacKay [Mac93j. 

An orbit ( |19[ ) is symmetric with respect to a reversor S if S(T) = T; i.e., there exists a j G Z 
such that S(xo, Zq) = (xj, Zj). Denoting the fixed sets of the reversors by Fix(Sj), it is well-known 
that any symmetric orbit necessarily has points on these sets, see App. [Aj Indeed, for the lift, 
there are two points on two fixed sets that occur essentially halfway around the orbit (depending 
on whether n is even or odd), and essentially halfway along in the angle direction (depending on 
whether the components of m are even or odd) , see Table [l] For the reversors (|25|) and translation 
symmetry (20), the relevant fixed sets are the curves 



Fix(Si o T- m ) 
Fix(S , 2 o T_ m ) 



{( 



2 ' - 



z 6 M}, 



{a^l{z,5) + m),z) 



: z G 



(26) 



Projecting back to the torus T 2 , we see there are eight symmetry curves — four for Si, the lines 
above x = (0, 0), (h, 0), (0, |) and |), and four curves for 5*2. Since symmetric orbits intersect 
these lines in pairs, as shown in Table [TJ we expect to find four distinct symmetric periodic orbits 
for each (m, n). 

Moreover, since the one-dimensional fixed sets are graphs over the action coordinate, z, any 
rotational invariant torus of / necessarily intersects all of these fixed sets. Consequently, there 
are symmetric orbits on each rotational torus, and if the dynamics on the torus is conjugate to an 
incommensurate rotation, then the symmetric orbits are dense. Therefore using symmetric periodic 
orbits as limiting approximations to rotational tori seems reasonable. 
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3.3 Stability 



The linear stability of a period-n orbit, (x*, z*) = f n (x* , z*), is determined by its Jacobian matrix, 
Df n (x* , z*). When / : M — >• M is volume-preserving and M is three-dimensional, the characteristic 
polynomial takes the form 

det(A/ - Df n (x*, z*)) = A 3 - riA 2 + r 2 A - 1, 

where the trace T\ and second trace r 2 are defined by 

ri=tr(Df), 
r 2 = |(r 1 2 -tr( J D/ 2 ")). 

There are eight generic multiplier configurations in the ti-t 2 plane as shown in Fig. [2] [L M98] . The 
line t\ = r 2 corresponds to the existence of a multiplier A = 1. The only (ti,t 2 ) values corresponding 
to stable orbits lie on the segment — 1 < t\ = r 2 < 3 where there is a conjugate pair of multipliers 
on the unit circle. When t\ > 3, this pair becomes hyperbolic, and, when t\ < —1, reflection 
hyperbolic. 



(27) 



1<X,=X, 




Figure 2: Stability diagram for a three-dimensional, volume-preserving map depending upon the trace Ti and second 
trace T2. The eight insets show multiplier configurations in the complex- A plane relative to the unit circle for each of 
the eight stability domains. 



The Jacobian Df of ( 16 ) is the matrix (in 2x1 block form) 



Df 



I - eVn(z')Dg(x) VO(z r 



-eDg(x) 



1 



(28) 



In this expression, Dg is the row vector (Vg) T , while VJl(z) is a column vector, so that V£lDg is 
the outer product of these vectors. 
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Interestingly, (28) has the property that, for n = 1, 



n = r 2 = 3 - £VQ(V) • Vg, 



(29) 



for any functions and g. Consequently one of the multipliers of Df is always 1, and thus every 
fixed point has multipliers (1, A, 1/A). For example, in Fig.[T]the unstable fixed points are indicated 
by dashed (t± < —1) and dotted (ri > 3) curves. For the parameters (18), the fixed points on the 
line x\ — X2 = \ are stable up to e = [tt(2(3z* — with z* given by (22). At this parameter 

value, ri = — 1 for the point (|,0, z*), and as e grows there is an interval of unstable fixed points 
orbits on this line. Similarly on the line X2 = the point (0, 0, z*) first loses stability with t\ = —1 



at e 



The fixed points on the line X2 = q have t\ > 3 for all e > 0. 



3.4 The Residue 



We noted above that n = r 2 for fixed points of any map of the form (16). Of course the equality 
r\ = r 2 need not be true for Df n (x, z), so periodic orbits with periods larger than 1 may not have 
a unit multiplier. However, for symmetric orbits, the property n = r 2 follows from the conjugacy 
of Df to Df' 1 . 

Lemma 2. IfT is a symmetric periodic orbit of a reversible, three-dimensional, volume-preserving 
map, then it must have a multiplier A = 1, or equivalently, t\ = r 2 . 

Proof. This follows from the more general result: each symmetric periodic orbit of a reversible 
diffeomorphism has reciprocal multipliers, i.e., if A is multiplier, then so is A -1 . Indeed by Lem. [4] 
in the Appendix, there is a point x* £ T and a reversor S such that x* £ Fix(S'). Moreover, since 
S*o/ = f'l-oS, then f n = S^of-^oS, so that the matrix A = Df n (x*) is conjugate to Df~ n {x*). 
Since Df~ n (x*) = A -1 , A and A~ l have the same spectrum. 

Finally, for the three-dimensional, volume-preserving case, the product of the three multipliers 
is one, and since the multipliers come in reciprocal pairs, there must be a multiplier A = 1. This is 
equivalent to T\ = r 2 . □ 

Using this result, we let r = t\ = t 2 and define an analogue to Greene's residue ^ for symmetric 
orbits 

R=\(3-r), T = tr(Df n (x,z)). (30) 

This definition was also used for a study of periodic orbits of the rotating standard map [ACS92, 
ITom96b] — reversibility was not required for that map, since one of the three multipliers is identically 
equal to one for all orbits. 

Just as for Greene's residue, orbits are stable when < R < 1 where they have a complex 
conjugate multiplier pair on the unit circle; these are analogous to elliptic orbits. The complex pair 
hits minus one at R = 1 and is subsequently replaced by reciprocal pair of negative multipliers when 
R > 1. When R = the orbit has a triple multiplier A = 1, and when R < there is a reciprocal 
pair of positive eigenvalues; these orbits are analogous to hyperbolic orbits of an area-preserving 
map. 

Analytical formulae for the positions and residues of the symmetric fixed points are shown in 
Table [2} The fixed points lie on Fix(Si) n Fix(5 2 ), recall ([26]), since m = £l(z*). Their residues 
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depend upon the parity of the components of m, which we denote by "e" for even and "o" for odd. 
There are four parities as shown in the table. 

Symmetric, period-two orbits must solve the system (24). Thus on Fix(Si), the orbit is 

mi 



(0. 



,m mi 



7), 



for 5* = \/3(m\ — 2j) 2 — Since gcd(mi,m2,n) = 1, there are three types of such orbits, 
depending upon the parity of m. Their residues depend similarly on this parity, see Table [2] A 
period- two orbit that begins on Fix(S < 2) is given by 

(i«(z*),z*) i — y (m-±n(z*),mi-2rf-z*), 

where the initial action must satisfy the transcendental equation 

2z* - = mi - 2 7 , 



and 5* is determined by (24). The residue of these orbits does not have a simple analytical form. 



(m, n) 


(or*,**) 




5* 


R 


(e,e, 1) 
(e,o, 1) 
(o,e, 1) 
(o,o, 1) 


(0,0, mi - 
(0, |,mi - 
(|,0, mi - 

( 5 ,5' m i " 


7) 
7) 
7) 
7) 


pz* 2 - m 2 


7T£ 



7re(2/3z* - 1) 
-27r/3ez* 


(o,o,2) 
(o,e,2) 
(e,o,2) 


(0,0,^- 
(0,0,^- 
(0,0,^- 


7) 
7) 
7) 


(3z* 2 - \m 2 


2vre(l - 2pz{l - 2vre)) 
27re(27re + 2/3z(l -2?re)) 
2vre 



Table 2: Properties of the four types of fixed points and three types of period- two orbits beginning on Fix(Sa) for 
a = b — c = 1. Orbits are classified by the parity of m. 



4 Computing Symmetric Orbits 

The computation of symmetric periodic orbits can be reduced to finding the zeros of a function 
H : M 2 — > M 2 since the orbits must begin and end on fixed sets of the reversors. Indeed, for a given 
fixed set, the initial angles are completely determined by (z,5). We therefore need only determine 
zq = z* and 5* by requiring that x n = xq + m. In addition, the numerical work can be halved: 
since the "half-orbit" lands on another fixed set, recall Table [TJ only \ ~ n /2 iterations need to be 
performed. 

More concretely, let (x l (z,S), z) £ Fix(S'j) and {x* (z, S), z) 6 Fix(Sj) denote points on the 
initial and final symmetry curves of an (m, n)-orbit, given in Table [l] using (p6|). If, for example, 



Sf = S\o T_ m then x* (z, 5) = |m, and if Sf = S 2 o T_ m then (z, 5) = ^(^(z, 5) + m). Denoting 
the length of the half-orbit by x so that (x x , z x ) = F x (x' l (z, 5), z), then an (in, n)-orbit corresponds 
to a zero of 

H(z,5)=x x (z,6)-xf(z x (z,6)), (31) 
that is, to a point (z*,S*) such that H(z*,5*) = 0. 
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Numerical solution of this two-dimensional system is straightforward using a Newton or pseudo- 
Newton method; we employ Broyden's method [DS96j . Recall that finite difference approximations, 
such as those used in the secant method, are underdetermined in more than one dimension. The idea 
of Broyden's method is to begin with an initial approximation to the Jacobian, the identity in our 
case, and improve it at each iteration by taking the solution of the finite difference approximation 
that minimally modifies (in the Frobenius norm) the current approximation. Using double precision 
arithmetic, we generally define convergence as \H(z*,6*)\ < 1(T 10 ; however, orbits were often found 
to higher precision. A maximum of 75 Broyden iterations were performed. 

An example of the dependence of (z* , 5*) on e for three low-period orbits is shown in Fig. [3j note 
that the dependence on e appears to be smooth. Orbits with periods up to O(10 5 ) and moderate 
values of e can be easily found. Moreover, we are able to find symmetric orbits for any coprime 
(m,n) and any pair Si, Sf as listed in Table [TJ 

It is important to note that errors can be compounded in these calculations. While the half-orbit 
error was always smaller than 10 -10 , the error for the full orbit, 

« = \%n - xq - m\ + \z n - z \, (32) 

was often larger, but generally no more than 10~ 8 . Of course when the orbit was stable or nearly 
stable, R = 0(1), the error k was typically smaller. 

Unstable orbits are difficult to find because the derivative of H becomes large and the algorithm 
becomes extremely sensitive to the initial guess. This can be partially obviated by using continua- 



tion from the trivial case e = where the initial conditions are given by ( 22 ) . We found quadratic 
extrapolation to be sufficient to predict subsequent initial conditions. Extrapolation typically re- 
quires 10-15 fewer iterates of the Broyden method, allows larger steps in e, and converges for larger 
residue values. 

To maintain convergence, the extrapolation step size must decrease as the orbit period grows. 
To start, the simple choice As = O.ln -1 yields few convergence failures. If successful, the step size 
is increased by half. If unsuccessful, the step size is decreased by a third. Though this process 
generally results in more failed steps than successful ones, it proved far faster than a method with 
a fixed step size. 



5 Computing Residues 

Given numerical values of (z*,5*) for an (m, n)-orbit, the residue ( |30[ ) is obtained by multiplying 
the successive Jacobian matrices along the orbit to compute ti(Df n ). Every (to, n)-orbit has zero 
residue at e = and we observe that R grows or decreases smoothly with e. The behavior is initially 
monotone, but for intermediate values it often oscillates, as can be seen for several low-period orbits 
with (x l ,z l ) G Fix(Si) in Fig. |4j For example, the residue of the (2,1,4) orbit (red curve) initially 
becomes negative, reaching a minimum near e = 0.04505 and subsequently grows monotonically, 
passing R = 1 near e = 0.1146. The residue of the (1,3,4) orbit (black curve), on the other hand, 
is initially positive, but reaches a maximum near e = 0.09815, and then decreases monotonically 
thereafter. The residues of higher period orbits can oscillate several times, but always appear to 
grow monotonically for large e. 

Note that the residue for fixed points and period-two orbits either grows linearly or quadratically 
for small e (except for the (e, o, 1) orbits which have zero residue), recall Table [2j Surprisingly, we 
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0.5 




Figure 3: Numerically computed (z* , 8*) for three low-period orbits beginning on Fix(Si) as a function of e using 
a step size of 0.005 




Figure 4: Residue as a function of e for four (m, n)-orbits of (16 1 that start on Fix(Si). The bifurcation values 
R — and R = 1 are indicated by dashed lines. 



observe that linear or quadratic growth holds for higher periods as well, some samples are shown 
in Fig. [5] There does not appear to be a simple relation between the period or type of orbit and 
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the small-e asymptotic behavior. For example, we observed the growth rates: 

f 0(e) : (1,1, 3), (2, 2, 3), (2, 3, 3), (3, 2, 3), (3, 3, 4), (1,4, 4), (3, 4, 4), (4, 3, 4), (1,1, 5),... 
\ 0{e 2 ) : (1,2, 3), (2, 1,3), (2, 1,4), (1,3, 4), (3, 2, 4), (1,3, 5), (3, 2, 5), (4, 2, 5), (4, 3, 5),... ' 

when £<1, independent of the initial symmetry lines. 

This asymptotic behavior contrasts with the behavior of the standard map, where R = 0(e n ) 
for e < 1 |Gre79] . An explanation for our observations could follow from a generalization of the 
Poincare-Birkhoff theorem to the volume-preserving case, but we know of no such result rj 




Figure 5: Residue as a function of e for three orbits of (161 that start on Fix(Si) plotted on a log-log scale. The 



cusps correspond to points where R changes sign. The asymptotic behavior is indicated by the line segments. 

Greene also observed that for the standard map, R = 0(e n ) for e > 1, We similarly observe 
that most of the residues of (16) grow as e n , three examples are shown in Fig. [5j However, there 



are exceptions: every symmetry line has a class of orbits — those listed in Table [3] — for which the 
residue grows as e n ~ l . Moreover, on Fix(Si o Tr^i) the residues of all orbits have this lower growth 
rate. 

As we observe in Fig. [4j the sign of the residue varies with the orbit and with e. Greene and 
MacKay found that the standard map has a special symmetry fixed set, the so-called dominant 
set, on which every orbit of the standard map has positive residue [Mac92j. This holds as well 
in other two-dimensional examples |RQ92| . To our knowledge no such result is known for higher- 



1 For example, Cheng and Sun's generalization of Poincare-Birkhoff applies to invariant circles, not periodic orbits 
|CS90b] , 
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Si 


Si o T 0) i 


Si o Ti i0 


Si T 1>t 


s 2 


S2 To,! 


S2 Ti j0 


S2 Ti ( i 


(e,o, e) 


All 


(0,0, e) 


(o,e,e) 


(e,o, 0) 


(e,e,o) 


(0, 0, 0) 


(0, e, 0) 



Table 3: Classes of orbits for which R = 0(e n x ) for £ > 1, grouped by symmetry line. Orbits not listed have 
R = 0{e n ) for e > 1. 



dimensional maps, though there is some numerical evidence of a dominant symmetry for the 4D 
Froeshle map |KM89j . 



To determine whether the map ( 16 ) has a dominant symmetry, we looked at a number of (m, n)- 
orbits of each symmetry and parity type. The results are summarized in Table |4j which indicates 
the stability of the orbit on each of the eight symmetry fixed sets. Unfortunately, we were unable 
to find any systematic organization of residue signs on any fixed set. One pattern that seems to 
hold is that for small, positive e, there appear to be two elliptic, < R < 1, and two hyperbolic, 
R < 0, orbits for each (m, n); these are labeled and hi for i = 1,2 in Table |4j This does not 
hold, however, for large e: the asymptotic sign of the residue as e — > 00 is indicated by the ± signs 
in the table. 



Orbit 


Si 


Si r 0> i 


Si Tifl 


Si Ti,! 


s 2 


s 2 r 0i i 


S2 Ti )0 


S 2 Ti,! 


(2,2,3) 


ei- 


hi- 


h 2 - 


e2+ 


ei- 


^1- 


^2- 


e 2 + 


(4,2,5) 


hi- 


h 2 + 


e 2 + 


ei- 


/11- 


h 2 + 


e 2 + 


ei- 


(2,1,4) 


hi+ 


hi + 


ei+ 


ei + 


/12- 


h 2 - 


e 2 - 


e 2 - 


(4,3,4) 


ei+ 


ei+ 


h 2 - 


/i2- 


h x + 


hi+ 


e 2 - 


e 2 - 


(2,1,3) 


ei+ 


hi+ 


e 2 - 




h x + 


ei+ 


h 2 - 


e2- 


(4,3,5) 


ei- 


hi+ 


e 2 - 


h 2 + 


h x + 


ei- 


h 2 + 


e 2 - 


(1,4,4) 


h x + 


ei+ 


hi+ 


ei+ 


e 2 - 


/* 2 - 


e 2 - 


^2- 


(3,2,4) 


ei + 


hi+ 


ei+ 


/ii+ 


h 2 + 


e 2 + 


h 2 + 


e 2 + 


(1,2,3) 


hi+ 


h 2 - 


ei+ 


e2+ 


ei+ 


e 2 + 


hi+ 


/l 2 - 


(3,2,5) 


ei+ 


hi + 


/i 2 + 


e 2 + 


h 2 + 


e 2 + 


ei + 


/ll + 


(1,3,4) 


ei- 


e2- 


e 2 - 


ei- 


h x - 


/i 2 - 


h 2 - 


/ll- 


(3,3,4) 


ei- 


hi- 


hi- 


ei- 


e 2 + 


/12- 


h 2 - 


e 2 + 


(1,1,3) 


ei + 


hi- 


h 2 - 


e2- 


e 2 - 


^2- 


h x - 


ex + 


(1,3,5) 


hx+ 


e2- 


ei- 


e2- 


e 2 - 


ei- 


e 2 - 


hx+ 



Table 4: Signs of the residues of symmetric orbits of each parity on each of the eight symmetry fixed sets for ( 16 1 
with parameters (18 1. Orbits are labeled by stability for small positive e ("e" for elliptic and u h" for hyperbolic), 
and by the residue sign (±) for large positive e. For small e there are two hyperbolic and two elliptic orbits for each 
(m,n), labelled by subscripts 1 and 2. 



The accuracy of the residue computations is less than that of the orbit itself because multipli- 
cation of the Jacobian matrices along the orbit gives rise to additional error. Roughly speaking, 
for fixed matrices A and B and k <C 1, 

tr((A + nB) n ) « tr(A n ) + ©(n^™" 1 ) 

where /i is the largest eigenvalue of A. For the residue computation we can think of fi as an estimate 
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of the Lyapunov multiplier or mean residue of the orbit, and thus the relative error in the residue is 
of order If k were fixed and \x = 0(1), the error would grow linearly with the period. However, 
as the period grows, the orbit error k, (32), increases because of the floating point errors in the 



computation of (31). 



This error estimate for the residue computations is shown in Fig. [6^a) for a period 5842 or- 
bit. When e < 0.023, \R\ < O(10~ 7 ), and the error makes the residue essentially uncomputable; 
however, the main point is that all of the multipliers for small e are essentially 1.0. When the 
orbit loses stability, near e = 0.025, it does so rapidly reaching R = 10, 010 at e = 0.0265, and the 
estimate indicates that when 10~ 7 < \R\ < 10 4 , the relative error in R is small. Note that even 
when R ~ 10 4 , the Lyapunov multiplier of this orbit is [i ~ R x l n = 0(1); consequently the error in 
the residue is still of order nn. 

The residue as a function of period for 35 orbits for e = 0.01 is shown in Fig. [6^b), together 
with the error estimate. This particular sequence of orbits corresponds to the Farey sequence of 
approximants of the "spiral mean" rotation vector that will be studied in §6.3| For this e, these 
orbits become increasingly stable as the period grows. While we are not able to obtain an accurate 
value for the residue when n > 300, we can nevertheless be sure that it is very small: the multipliers 
of these orbits sit very close to the point r = 3.0. We will argue in {7^2 that this sequence of periodic 
orbits appears to converge to an invariant torus. 



(a) 



(b) 



i - 



\R\ 
hi s 



10- 5 ■ 
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Figure 6: (a) Residue of the (1897,4410,5842), Fix(Si) orbit as a function of e. (b) Residues of a sequence of 
Fix(Si) orbits with periods 1 to 17991 (the orbits of Table [bJ with e = 0.01. The bars indicate the estimated error 



6 Farey Sequences 

As we discussed in ^2j the application of Greene's residue criterion requires finding a suitable 
sequence of (m, n)-orbits whose rotation numbers converge to a chosen Diophantine rotation vector. 
For the d = 1 case, these can be systematically obtained either using continued fractions — for which 
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the successive convergents Q are "best approximants" — or the Farey tree — which gives, in addition 
to the convergents, intermediate approximants. 

There is no multi-dimensional generalization of the continued fraction algorithm that is guaran- 
teed to give exactly the sequence of best rational approximants to a given irrational vector, though 
there are strongly convergent algorithms |Lag94, KDM07]. In lieu of these, we will use Kim and 
Ostlund's generalization of the Farey tree to generate rational approximations to incommensurate 
rotation vectors [K086 . 

In addition, it is not at all clear what Diophantine vector may be an appropriate choice to replace 
the golden mean, (10), which appears to give the most robust invariant circle in the standard map. 
Since, as we recall below, integral bases of degree-(d + 1) algebraic fields can be chosen to make 
Diophantine vectors, it seems natural to conjecture that some cubic field will lead to robust two-tori. 
We will investigate vectors from the so-called "spiral mean" field. 



6.1 Generalized Farey Tree 

The Farey tree is a recursive algorithm for generating all rationals between a given initial pair 
[HW79J. To do this, begin by regarding each uj G M + as the direction (uj, 1) t G and each 
positive rational ^ as the direction (m,n) T . Conversely, for each vector (u,v) T G there is an 
associated rotation number uj = u/v. The problem of approximation is to then find a sequence of 
vectors (m£,n£) T that approach (uj, 1) t in the sense that the two directions become parallel, i.e., 
\ri£UJ — mil ~ * as i — > oo. 

The root of the Farey tree consists of the two vectors, (0, 1) T and (1, 0) T , said to be at level zero. 
The interior of the cone of directions spanned by the positive linear combinations of these vectors 
thus represents all positive numbers. The set of directions can be mapped onto the segment of the 
line u + v = 1 in the positive quadrant, and each direction in the "level-zero" cone corresponds to a 
point on this segment as sketched in Fig. [7| The next level is obtained with the mediant operation 
©, 

( m )®( P ) = ( m l P ), 
\nj \qj \n + qj 

which is simply vector addition applied to the bounding vectors. Thus at level one there is one 
new daughter vector, (1, 1) T = (0, 1) T © (1, 0) T . The daughter divides the level-zero cone into two 
"level-one" cones 

C/= (i i)> a= (o l)' 

where the columns represent the bounding vectors. This process can be repeated: each level-^ cone 
is divided by forming the mediant of its boundary pair. Thus the level-one cone C\ becomes the 
two cones C/C; and CiC r , and there are 2 £ cones at level I. 

The Farey tree is then the binary graph obtained by connecting each rational vector at level 
I > 1 with its two daughters at the next level. Each vertex has a unique path, p G {/, r} N : the path 
to (1, 1) T at level one is empty, and p = I for the daughter (1, 2) T to the "left" and p = r for the 
daughter (2, 1) T to the "right". Each level-£ parent has two level- (I + 1) daughters, whose paths 
are obtained by appending an / or r, respectively, to the parent's path. The vector corresponding 
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to a path, p = i\ %i . . . ig, € {I, r} is then 

(n) = C ix C *-~ C it (J) 




/ \ r' V , r' \ / V 

(1,4) (2,5) (3,5) (3,4) (4,3) (5,3) (5,2) (4,1) 



Figure 7: First five levels of the Farey tree construction as the set of directions in the first quadrant projected onto 
the segment x + y = 1. 

Since the matrices C\ and C r are unimodular, the parallelogram formed from the two vectors 
bounding each cone has unit area. A consequence is that every vector on the Farey tree is coprime, 
gcd(m, n) = 1. Moreover, each rational has a finite Farey path, and every irrational to is obtained as 
the limit of an infinite path on the tree. For example, the golden mean has path <j) = rlrlrlrl . . . = rl. 
There are also infinite paths that limit on each rational from above and below, these are the paths 
with infinite tails of Z's and r's, respectively. There is also a simple relation between the continued 
fraction expansion and the Farey path, see e.g., [HW79J. 

Kim and Ostlund generalize the Farey tree to higher dimensions by regarding each vector 
lo £ Rl as a direction (to, 1) T £ JK086]. We will take the level -zero cone to be the positive 
octant, generated by positive linear combinations of the triplet (1,0, 0) T , (0,1, 0) T , (0,0, 1) T . The 
intersection of the cone (u,v,w) T G with the plane u + v + w = lis the equilateral triangle 
sketched in Fig. [8](a). Somewhat arbitrarily, the level-one child is defined to be the mediant of two 
of the bounding, parent vectors; we choose (1, 1, 0) T = (1, 0, 0) T ® (0, 1, 0) T . This creates a pair of 
level-one cones bounded by the child and the two parents that form an independent triplet, that is 
the cones 

/0 1\ /0 1 1\ 

Ci = 1 1 ) , C r = 1 . (33) 
\0 10/ \1 0/ 

As shown in Fig.^a), the line connecting the unused parent, (0, 0, 1) T to the child, (1, 1, 0) T divides 
the equilateral triangle into "left" and "right" triangles. To continue the process, a new child in 
each cone is obtained by applying the mediant operation to the two oldest parents (those generated 
from the earliest levels), and each cone is divided into two, using the resonance line connecting the 
child and its youngest parent. The construction up to level three is shown in Fig. [8^a). 

Projecting the cones onto the plane w = 1 gives the frequency plane, as shown in Fig. [8)3. Each 
point to in this plane represents a direction (to, 1); here we show only rotation vectors in the unit 
square [0, 1] x [0, 1]; these correspond to directions in the cones rrl and llr of the full tree. 
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Figure 8: (a) First three levels of the generalized Farey tree projected on the plane :r + y + z = lin homogeneous 
coordinates. The two shaded triangles correspond to rotation vectors in [0, l] 2 . (b) Projection of the two level-three 
cones, rrl and llr, onto the plane w = 1, and three additional levels of the tree. 



This binary, recursive procedure assigns a path p £ {I, r} N to each vertex below level one by 
following the sequence of dividing lines from (1, 1, 0) T : a turn to the "right" appends an r to the 
path and a turn to the left, an I. Thus, for example the path rrr corresponds to (2, 1, 1) T and Irl 
to (0, 2, 1) T . Just as for the Farey tree, the vector corresponding to a path p is 



(34) 





using the matrices (33). 



Since the matrices Cj are unimodular, the three vectors that bound each cone define a paral- 
lelepiped with volume one. A consequence is that every vector on the tree is coprime, gcd(mi, m2,n) = 
1. Moreover it can be shown that every coprime vector in the positive octant has a finite path on 
the tree. One defect of this construction, however, is that every integer vector in the interior of the 
positive octant has two finite paths; for example, rrlll and rrlrr both lead to (2, 1,3) T {^] Just as 
for the Farey tree, however, every incommensurate direction has a unique, infinite path. There are 
also infinite paths that limit on rational directions; for example, any path p = hlr with a periodic 
Ir tail and arbitrary head h. 

Accounting for the duplication, there are 

2 £ + [5-(-l)']2LU- 1 + l 



2 If (m, n) lies on a level-^ dividing line then its two paths are those obtained by exchanging I and r for all entries 
beyond level I. 
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Figure 9: (a) The 2140 unique rotation vectors that give (cji,^) £ [0, l] 2 on the Kim-Ostlund tree up to level 13. 



(b) The corresponding values (z*(u;), <5*(w)) for the (m,n)-orbits of (161 for e — 0. Rotation vectors on the four lines 
^1,2 = 0, 1 are colored and labeled in (b). 



unique directions on the tree up to level i (including the initial triplet). In Fig. [9^a) we show the 
2140 unique rotation vectors in the rrl and llr cones up to level 13. The boundaries of each cone 
at each level correspond to resonances p ■ lv = q, and these give rise to the isolated lines in the 
figure. Just as for the Farey tree, the sampling of rotation vectors is nonuniform, tending to avoid 
low-order resonances. 

Each of the rotation vectors in Fig. [9]^a) corresponds to a periodic orbit of (16) at e = 0. The 
(z*,5*) values, given by (22), for these orbits are shown in Fig. [9^b). Note that the transformation 
uj h- > (z*,5*) turns each line <jJ\ = const into a parabola, and that since the map £l(z,5), (12), is 
orientation reversing, the upper and lower boundaries of the square are flipped. 

Each of the points in Fig. [9^b) continues to an (to, n)-orbit of the map (16) on each symmetry 
line for c / 0. The orbits on Fix(Si) for e = 0.02 are shown in Fig. 10 The points are colored 
to indicate their residues: orbits with the smallest residues are light (yellow), those with \R\ ~ 1 
are dark (red), and those with \R\ > 1 are black (small points). The corresponding (z*,5*) values 
for each of these orbits, Fig. [lO^b), are shifted from their e = values so as to widen the gaps 
around the low-order resonance curves. In phase space this correlates with the growth of tubes 
of secondary tori for each resonance, and these are largest for the forced resonances, p = (0, 1), 
(1,0), and (1,-1), of the force (15) [D M12] . Similar behavior has been observed in action space 
for symmetric orbits of the four-dimensional Froeshle map [KM89 . 

The lightest regions in Fig. 10 correspond to groupings of the most stable orbits. As we will see 
m ^ these regions also contain the most robust tori. 
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Figure 10: Periodic orbits of (16 1 on Fix(Si) for e = 0.02 with rotation vectors as in Fig. [9] The color indicates the 
magnitude of the residue, \R\ < 1, as shown in the color bar. Orbits with \R\ > 1 are shown as small, black points. 



6.2 Cubic Fields 

KAM theory suggests that the most robust invariant tori should have Diophantine rotation vectors. 
One way to construct these is to use integral bases of algebraic fields. Recall that a G C is an 
algebraic number of degree n if it is a root of an n th -degree polynomial with rational coefficients, 
but not of any such polynomial of lower degree. The algebraic field generated by a is 

Q(a) = {Zj *"' : Oj 6qJ- 

For example, the golden mean c/>, the larger root of c/> 2 — c/> — 1 = 0, is a quadratic irrational and 
Q(t/>) = {a + b(j) : a, b G Q}. The ring of integers Z(a) C Q(a) consists of those elements that are 
roots of monic polynomials, x n + p\x n ~ l + . . . + p n , with pi G Z. The ring is an n-dimensional 
vector space over Z, and a basis {a±, . . . , a n } for this space is called an integral basis. In general, 
integral bases give rise to Diophantine rotation vectors. 

Theorem 3 ( |Cas57| ICus72j ) . // (w, 1) G ~R. d+1 is an integral basis of a real algebraic field of degree 
d+X, then uj G T>d- 

For example, (</>, 1) is a basis for Q(c/>)> thus Z(</>) = {q<p + p : q,p G Z}. Every integral basis 
for Z(c/>) is of the form (01,02) = (</>, X)M, where M G SX(2,Z). The ratios, u = 01/02 of these 
integral bases are precisely the noble numbers, conjectured to correspond to the locally most robust 
invariant circles of twist maps, recall §[T] 

Thus a natural extension of the robustness conjecture of the nobles for d = X is that the most 
robust tori for d = 2 should have rotation vectors in an integral basis of a cubic field. Several 
candidates have been proposed for the most robust field K()8(i. Lo_c92, Tom96aj, and invariant 
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tori with such rotation vectors have been studied for 4D symplectic maps, recall ^T} However, the 
robustness conjecture has not yet, to our knowledge, been investigated. Indeed, it would be difficult 
to find a "last" invariant torus in the symplectic case: since the tori are not codimension one, they 
do not form barriers. 



6.3 Spiral Mean 

One natural generalization of the golden field is that of the spiral mean, the real root of 



a 



a 



1 = 0, 



(35) 



proposed by Kim and Ostlund [K086J. We will denote the spiral mean by a and its conjugates by 
o\ and 0"2 = d\\ here 



a = a 
1 



1.3247179572447460, 



f 1,2 



a 



9 = arccos(-lai) rj 2.437734932288317. 



(36) 



The spiral mean, like (ft, is a "Pisot-Vijayaraghavan" (PV) number: a root of a monic polynomial 
with exactly one root outside the unit circle. Indeed, 4> and a are the smallest quadratic and cubic 
PV numbers, respectively. The discriminant of (35) is square-free, A = —23, so the set {a 2 , a, 1} 
is an integral basis for Z(<r). 

For our purposes, it is providential that Z(<r) is naturally related to the generalized Farey tree 
construction. Indeed the polynomial (35) is the characteristic polynomial for the cone matrices Ci 
and C r of (33), and their maximal eigenvectors, 



(37) 











H 




, Vf=\ 













are in Z(cr). Moreover, since every entry of the fifth powers of both matrices (33) is positive, 

t>3 



these matrices are irreducible in the Perron- Frobenius sense. As a consequence, for any v G 



Cfv — > ca e Vi as £ — > oo. Thus by (34), the directions vj and v? correspond to the paths p = l°° = I 
and r°° = f, respectively. 

As an example, the sequence of integer vectors (mi, ni) with paths p = r e are 



(38) 





It is not hard to see that the sequence of periods {no, n\, . . .} satisfies the recurrence 

n i+3 = n t+ i + n t , n = 0, n\ = n 2 = 1, 

which has the solution 

(.SI 
rig = a^a + a\a l + a^o^, 

a, = 2^(3(7? + 7(Ji - 2). 



(39) 
(40) 
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The directions of Farey sequence ( |38| ), mi/ rig = (n^+i/n^, n^_i/ra^), limit on the rotation vector 

ijj = (a,l/a) = (a, a 2 — 1). 

Unlike the equivalent sequence for the golden mean, these successive approximants include rotation 
vectors that are not "best" approximants to u. Recall that a rational m/n is a best approximant 
to cj if 

\\nui — m\\ = min min \\quj — p\\ (41) 

0<q<np£Z<i 

for some norm | ||. The sequence of best approximants depends upon the choice of norm and, 
unlike the classical Farey tree, it is not known if the generalized Farey tree gives all of the best 
approximants. However, up to period 50, 000, the best approximants to Vf in the 1, 2 and oo norms 
are the same; they have the periods 

n e e {1, 3, 9, 12, 37, 114, 151, 465, 1432, 1897, 5482, 17991, 23833}. 

Moreover, we found that all of the best approximants for these three norms can be found in the 



Farey sequence (38) up to at least period 10 6 . 



7 Residue Criterion for a Spiral Mean Torus 

An invariant torus on which the dynamics is diffeomorphic to rigid rotation must have each Lya- 
punov multiplier corresponding to its tangent directions equal to one. Consequently, all of the 
multipliers of a codimension-one torus in a volume-preserving map must be one. Thus a natu- 
ral conjecture is that the multipliers of a sequence of periodic orbits that limit on such a torus 
should all limit to one. Similar considerations for the symplectic case led to the residue bounds 
discussed in ^2] As a consequence, for the three-dimensional volume-preserving case, periodic orbits 
in the neighborhood of a smooth invariant torus should have stability parameters near the point 
r\ = T2 = 3 that corresponds to a triple-one multiplier, recall Fig. [2j 



When the map is reversible and the reversor fixed sets are graphs over the action (as in (26)), 
then every rotational invariant torus will contain symmetric orbits. If the torus has an incommen- 
surate rotation vector, then these orbits will be dense. Thus it seems sensible to investigate the 



existence of such tori by looking for sequences of symmetric periodic orbits whose residues, (30), 
approach zero. 

7.1 A Spiral Torus 

In this section we explore the generalization of Greene's residue criterion for one particular spiral 
mean torus. Here we study the direction 

llf ~ (1,<7 3 ,(T 4 ) T . (42) 

Projecting this direction onto the plane (w, 1) gives the rotation vector 

u = {a -I, a 2 -I) £ [0,1] x [0,1]. (43) 
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path 


(m, n) 


path 


(m, n) 


11 


(1,1,1) 


llr 2() 


(114,265,351) 


llr 


(0,1,2) 


llr 21 


(151,351,465) 


llr 2 


(1,2,2) 


llr 22 


(200,465,616) 


llr 6 


(2,5,7) 


Ur 33 


(4410, 10252, 13581) 


llr 7 


(3,7,9) 


llr 34 


(5842,13581,17991) 


llr 8 


(4,9,12) 


llr 35 


(7739,17991,23833) 


llr 9 


(5,12,16) 


llr 36 


(10252,23833,31572) 



Table 5: Integer vectors for the generalized Farey paths llr 1 



The Farey convergents for the path ( 42 ) are 



llr 1 ~ (n e ^ l7 ni + 2,ni +3 ) J 



where ri£ obeys ( 39 ) . Several of the vectors in this sequence are shown in Table [5] 
When e = 0, the llr 1 periodic orbits are, by (22), at 



n t-\ r o 2 n t+2 

7, S t = pz £ 



ne+3 ne+3 
These values converge geometrically to the position of the spiral torus at e = 

-0.2933160310, 



(44) 



Zoo 
Soo 



a -I 

■2 



7 



a' 



-0.5828090783, 



as can be seen in Fig. [lira). Indeed, from M0j n e ~ a a E + 0{a-^ 2 ), so the convergents approach 



this point geometrically at a rate a 3//2 ; for example, 

+ 0.4892) + 0(<j- 31 ) 



0.6037 

?00 H WTTTT^ COS( 



31/2 



(45) 



with 9 from (36). The line in Fig. 11 



a) corresponds to ( 45 ) with the cosine replaced by 1 ; it is an 



effective upper bound to the values (44). The modulation of this geometrical convergence by the 



phase is illustrated by the (green) curve in the figure, which is obtained from (44) by treating i 
as a continuous variable. The near periodicities of the discrete levels can be explained by the fact 
that « -j? and even more closely 9 w ^ . 

It is interesting that the geometric convergence of zg appears to be maintained for e > — 
computations of the actions of orbits for the spiral convergents for four values of e are shown in 
Fig. lib). In this figure, data is not shown for periods larger than 10 3 for the two largest values of 



e (0.026 and 0.0264) because the orbit finding algorithm did not converge. Note that even though, 
as shown in the inset, the converged point (zoo, <^x>) varies with e, the geometric convergence rate is 
consistent with the e = value o~ 3 l 2 \ however, there is some indication that the convergence rate 
could be faster for the larger two e's. Unfortunately, the oscillations make it difficult to extract a 
precise rate. As we will see below, for e > 0.0259 the torus seems to be destroyed, and by analogy 
with the standard map [MMP84], a different convergence rate might be expected. 
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Figure 11: The initial action for the first 35 approximants to ZZr°° for (a) e = and (b) four cases with e > 0. 



The curve (green) in panel (a) is the exact result ( 44 1 upon extending the level L to a continuous variable. The line 



(black) in both panels is the upper bound of the asymptotic result ( 45 1 with slope — § on the log-log scale 



7.2 Critical Torus 



Recall from Fig. |6jb) that the residues of the level-£ orbits appear to converge to zero for fixed, 
small s, but from Fig. [5] that each grows geometrically with e. As can be seen in Fig. [l2|a), the 
residues of low period spiral mean approximants begin to grow near e ~ 0.03, though for periods 
larger than 23, the onset of rapid growth in R occurs near e = 0.26. For larger period spiral 
approximants, as shown on a log scale in Fig. 12 b), rapid growth in R begins near e ~ 0.0259. 
This evidence suggests that the llr°° torus does not exist for e > 0.0259. 
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Figure 12: Residue as a function of e for Fix(S*i) orbits with rotation vectors to llr for (a) six low-period orbits, 



and (b) six longer period orbits near e cr . The stars indicate the residues for the three values of e shown in Fig. 14 
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Table 6: Stability thresholds for the first 36 spiral mean approximants of llr°° on Fix(Si). Shown are the values 
ER th at which the residue first reaches three thresholds, Rth = 0.5, 0.9, and 1.5. The * denotes orbits that are best 
approximants to w in the sense of (|41[). 



To more systematically investigate the residue criterion, we choose a threshold residue, R t h, 
and compute the smallest positive e for which the level I orbit has \R\ = Rth] call this value £R th - 
Computationally, we continue each orbit in e, as discussed in £j4j until Rth is bracketed. Bisection 
is then used to find £R th as precisely as possible. The threshold e values are shown for the first 36 
levels, up to period 31572, in Table [6} These thresholds appear to converge to 

e cr w 0.02590 ±5(10)" 5 , (46) 

independently of Rth- Computations for the spiral mean orbits on the other symmetry lines gave 
the same result. 



The convergence of the thresholds to e CT with period is geometric, as shown in Fig. 13 for 
Rth = 0.9. A least squares fit to this data gives 



£e = £cr + (0.17 ± 0.06)n- u -™^. (47) 

It is perhaps interesting that the rate of convergence, 0.79, is close to, but significantly different 
from, a' 1 0.75488. 
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Figure 13: The e values for R t u = 0.9 for the first 36 approximants to llr°° using e cr = 0.02590. 



7.3 Supercritical Tori 



These results imply that the llr°° torus appears to be critical at e cr 
Fix(Si) periodic orbits at this value is 



0.02590. The limit of the 



j OO ) ^OO ) u oo J 



(0,-0.300443,-0.581962), 



which appears to lie on on a critical spiral mean torus. 

By analogy with the area-preserving case, we expect that when e > e cr , the llr^ sequence will 
still converge to a quasiperiodic orbit, but one that is not dense on a torus; indeed, for twist maps 



the supercritical orbits are dense on Cantor sets. Figure 14 gives a visualization of this for our 
map: it shows the phase portraits of three spiral approximants for three values of e near e cr . The 
residues for these nine cases were marked with a * in Fig. 12 b). When e < e cr (the left column 



of Fig. 14 ) the approximating orbits are stable and as the period is increased the points appear to 



limit to a smooth, nonzero density on a smooth torus. In the middle column of the figure, where 
e = e cr , the density becomes less smooth, and in the right column, where e > e cr , it appears to 



limit to a density with holes. This can be better seen in Fig. 15 which shows the projection of the 
I = 34 approximant onto the (x, y) plane. 



7.4 Varying Parameters 



Up to this point we have studied only the standard parameter set (18); in this section we vary 



the parameter c in the force (15). 

(x,e,c) — > (x + (jj 2)5 ~~ e > ~ c )- Since we have observed that e 



Note that the map (16) is invariant under the replacement 

is the same for the Fix(Si) and 



the Fix(Si o T\ t \) orbits, the critical graph (e cn c cr ) will be symmetric under reflection through the 



origin. Thus we study only positive values of e. 

, C r 



Fi gure |16| shows the critical set {e c 



of the spiral torus (42) for c G [—2,2], holding the 



remaining parameters fixed according to (18). On this scale, e cr appears to be a graph over c, 
and for c < this graph appears to grow smoothly as c increases. The peak value, e cr w 0.06156, 
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Figure 14: Phase portraits of the periodic orbits ZZr 30 , llr 32 , and llr 34 (rows) on Fix(5*i) for three e values (columns). 
For s < Ecr ~ 0.0259 the orbits appear to limit on a smooth torus. At e cr the density is nonuniform, and when 
e > e cr the orbits seem to limit to a torus with holes. 



occurs at c cr ~ 0.07281. Though the graph appears to decrease monotonically for larger c values, 
this is not the case: there are a number of cusps. These appear to be associated with an increasing 
number of oscillations in the function R(e) before it settles into its asymptotic growth, recall £[5j 
An enlargement near the most prominent cusp is shown in Fig. 16 b). 

The oscillations in R(e) make it harder to obtain accurate values for the critical parameters near 



each cusp. To resolve the cusp shown in Fig. 16 b), we used a longer period orbit (period 31572), 
and a larger threshold R t h = 10. Note that what appeared to be a single peak near c = 0.073 
becomes, in the enlargement, a number of cusps and several possible discontinuities. Each cusp 
corresponds to a change of the large e asymptotic behavior of R, and near the cusp the number of 
oscillations in the graph R(e) changes. We are not able to resolve the apparent discontinuities in 
Fig. 16 b) — it is possible that the critical set is not a graph over c in these regions. Several less 



prominent cusps are also visible in Fig. [T6|a); the local behavior near these is similar to the one 
shown in the enlargement. 
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Fi gure 15; Projection of the Fix(S'i) spiral approximant Z/r 34 , with (m, n) = (5842, 13581, 17991), onto the angle- 
plane for e — 0.0265. 




o.os 



Figure 16: Critical set (e Gr , c cr ) for the llr°° torus with the remaining parameters fixed as usual at ( 18 1, (a) Critical 
set computed using the llr 30 , period 5482, orbit with R t h = 0.9. (b) Enlargement near the peak e cr using the llr 36 
orbit with Rth = 10. 



8 The Last Torus 

The critical function, e cr (uj) is the value of e at which the torus with rotation vector oj is first 
destroyed upon perturbation from integrability, e = (here we assume that 5 = (z 2 ) — 0J2, recall 
(23)). When u is resonant, p ■ oj = q, the torus is typically destroyed for any e / 0; i.e., e cr (uj) = 
at resonance. However, KAM theory implies that e cr (^) > when u is Diophantine; thus, the 
critical function is highly singular. 

Our goal in this section is to compute e C r(w) and to identify the most robust torus of our 
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model (16). Ideally we would like to discover the analogue of the golden mean, which, as Greene 
conjectured, appears to give the most robust invariant circle of the standard map (g]). Though 
there is considerable numerical support for Greene's conjecture, it has never been proven]^] Strong 
numerical evidence also indicates that invariant circles with "noble" rotation numbers, 

uj = Q< f \ , with ad — be = ±1, 

cq> + a 

are locally robust for area-preserving maps [Per82, [MS92] . That is, in any interval of rotation 
numbers the most persistent invariant circle has a noble rotation number. Recall that each noble 
number is the projection of the integral basis (cuf) + b,c(p + d) of Z(<^>) onto (uj, 1). 

8.1 Critical Function 

Thus it is natural to conjecture that the most robust tori for a two-angle map may correspond 
to integral bases of some cubic field. As a preliminary investigation, we will compute the critical 



function for bases of Z(cr). Recall that Vf (37) is one such integral basis and that every other integral 



basis is obtained from Vf by multiplication with some element of SX(3,Z). To construct a sampling 



of these matrices we will use the 2 directions on the generalized Farey tree at level £, recall {6.1 



To each path h = i\i<i . . . it with ij G {r, 1} there is a corresponding matrix C^C^ . . . Ci e G SL(3, Z) 



using (33). Therefore, the direction 

Vp = Ci 1 Ci 2 . . . Ci e Vf , ij G {I, r}, 

is an integral basis for Z(cr) with path p = hf. We refer to h as the head and f as the tail of 
this path. As we remarked in §6.1| some of the finite paths at level £ on the generalized Farey tree 
correspond to duplicated directions; however, the 1 directions obtained by appending the infinite 
tail f are unique. 

The projection of each v p onto (uj p , 1) gives a rotation vector ui v whose components are Dio- 
phantine by Th. [3j Since the matrices C\ and C r do not generate SX(3,Z), the Farey tree does not 
generate all of the integral bases. Nevertheless, as was seen in Fig. [9j the resulting rotation vec- 
tors Up give a discrete sampling of frequency space that is sensibly more dense in the nonresonant 
regions. 

We estimate e cr (u p ) for each such invariant torus using the residue criterion for nearby periodic 
orbits. To give a reasonable approximation, but limit the computation time, we approximate the 
torus using the periodic orbit with path hr k , selecting the length, k, of the tail so that the orbit 
has period larger than 5000. Note that the k we choose depends on the head h. We then compute 
the value of e for which \R\ = Rth = 0.9 for each of the corresponding periodic orbits hr k . 

To sample the tori with oj G [0, l] 2 , we began by generating all of the level I = 13 heads that 
start with rrl or llr, namely those in the two triangles of Fig. [8^b). The result, shown in Fig. 17, is 
an approximation to e cr (uj) on the 2 x 2 10 = 2048 rotation vectors at this level. In this figure, the 
size of the dots (and their color, as shown in the color bar) is proportional to e cr (ujp). Note that 
there are several local peaks in the critical function, and that it tends to be small near low-order 
resonance lines. The most robust torus in this figure has the path p = r 2 / 6 r 3 / 2 r fc with k = 19; this 



3 It is known that the standard map has no rotational circles when 2ire > 63/64 |MP85j . 
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Figure 17: The critical function for (16 1 with the standard parameters (18 1. The values of e cr are low near the 



resonances, but grow into several peaks. The level-8 triangle that contains the most robust torus, r I , is shown 
(dashed) . 



(3901, 1718, 5268), or rotation vector uj p « (0.7405, 0.3261). For this torus, 



corresponds to (m, n) 
£ cr (ojp) = 0.0424. 

Hypothesizing that the most robust torus is nearby, we zoom-in by focusing on the level-8 
triangle 

f 2 
>3 

that contains this vector, i.e., the dashed triangle in Fig. 17 We now generate the 2 9 , level- 17 paths 
within this triangle, i.e., each path of length 17 that begins with r 2 ? 6 . Appending r k to each of these 
to give a periodic orbit with period larger than 5000, gives the zoomed-in view of the critical function 
shown in Fig. 18 a). The most robust torus of the 512 sampled vectors has e 
r 2 l 6 (rlr) 2 lr 2 r k with the k 




(oo p ) = 0.0479. It 

has the path p = r l P{rlr)Hr A r k with the k = 17 approximation (m, n) = (4276,2081,5806). We 
again zoom-in near the most robust torus, incrementing the level of the base triangle by three, and 
focusing now on the level- 11 triangle 



r I rlr ~ 




32 



The critical function for another 2 9 rotation vectors using this level- 11 triangle as the base is shown 
in Fig. 18 b). The most robust rotation vector at this stage has the path p = r 2 l®(rlr) 2 l 5 rr k with its 



k = 14 approximation (m,n) = (4670,2278,6353) and e cr (co p ) 
in the level- 14 triangle 



r 2 r(rlrY 
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6 3 7 
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which is outlined in Fig. 18 b). 
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Figure 18: Enlargements of the critical function for two "most robust" triangles, showing e cr for 512 rotation 
vectors in the spiral field that have heads nine levels down from their respective triangles, (a) Enlargement of the 



triangle r I (dashed) from Fig. 17 The level-11 triangle, r I rlr containing the most robust torus is also shown 



(solid), (b) Enlargement of the level-11 triangle from (a), and the level-14 triangle (small dashes) containing the most 
robust torus at this stage. The location of the most robust rotation vector at level-38 is denoted by oj* . 

Continuing in this fashion, we increment the level of the base triangle by three and then find the 
most robust torus among 2 9 rotation vectors within this triangle. The peak of the critical function 
increases with each zoom, leading to the level-29 triangle 

r 2 l 6 (rlr) 2 r(rl) 4 r 3 l 3 ~ I 85 119 163 I . (48) 

\234 328 449/ 

The most robust torus in this triangle has path p = r 2 l G (rlr) 2 r(rl) r 3 l 3 r 2 (lr) 2 r , giving (to, n) = 
(4800,2371,6531) so that u « (0.7350,0.3630). This orbit has initial condition 

(x*,z*,5*) ^ (0,0.123605791956645,-0.331972596409587). 

The peak critical value has converged to e max = 0.0522 ± 0.0005, to three significant figures. Note 
that the given error bound signifies the observed variation in the maximal e values for period 5000 
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orbits. Since we fix the period, this e max does not necessarily represent the true maximal value of e 
for infinite period. Nevertheless, this orbit is our best estimate, to period 5000, of the most robust 



torus for (16) for the standard parameters (18). 



To improve this result, we re-examined the level-29 triangle (48). As before, we use the 2 y level 
29 + 9 periodic orbits within this triangle, but now add a tail r k so that the period is at least 
25, 000. Of these 512 rotation numbers the most robust torus has the path 



p = r 2 l 6 {rlr) 2 r{rl) 4 r 3 l 3 r 3 (lrfr 8 , (49) 

which agrees for the first 31 levels with the previous one, but differs after that. This corresponds 
to (m, n) = (23749, 11731, 32316), so that 

u) rj (0.73490,0.36301). 

This torus is located at 

(x*,z*,5*) « (0,0.123303207885153,-0.332181389896200) 
and is destroyed at the critical value. 

e m ax = 0.0512 ± 0.0005. (50) 

This value of e m ax is less than the value of e max computed using the period 5000 orbits, which is 
consistent with the behavior seen in Table [6j Namely, critical e values, by and large, decrease as 
the period of the approximating periodic orbits increases, limiting on the true e cr from above. 
Unfortunately, unlike the noble numbers of the area-preserving case, there does not appear to 



be any simple pattern in the Farey sequence (49) of the most robust torus. 



8.2 Crossing Time 

To compare these estimates for the peak values of e cr with the dynamics we performed a "crossing 
time" experiment, similar to that of [Meil2j . The point is that a rotational torus is a barrier so 
that any orbit that begins below the torus must remain below. 

To make this computation easier, we choose a frequency map that is periodic mod-one in z, 
i.e., so that Q(z + 1) = Q(z) + m for some m 6 Z 2 . In this case, the invariant sets of the map 
in each unit interval of z are identical. Moreover there are no rotational invariant tori whenever 
there exists an orbit for which \zt — zq\ > 1 + A, where A is the maximal vertical extent of a torus. 



Since the first component of (12) satisfies the periodicity requirement, we modify only its second 
component: 

n 2 (z) = /3(nint(z)) 2 - S, 

where "nint" is the nearest integer function. Since the tori for parameters near (e, 5) = (0.05, —0.3) 
are located near z = 0.1 and have vertical extent of order 0.1, this modification will not affect their 
existence or nonexistence. Of course, it certainly changes the dynamics whenever \z\ > 0.5. 

The crossing time, t c , is defined to be the first time for which \zt c (xQ, zq) — zq\ > 1.3. If there is 
an orbit with t c < oo, the map has no rotational tori with A < 0.3. The crossing time is certainly 
a highly variable function of initial conditions, so we compute its distribution for a set of initial 
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conditions. Fixing 5 = —0.33, we start ten orbits with xq = and a random value \zq\ < 0.001. 
When e > 0.1 all of the orbits cross within a few thousand iterates, but as e decreases, the crossing 



times grow rapidly, as shown in Fig. 19 We found that whenever e < 0.0512 none of the trajectories 



crossed within 10 10 iterations. The implication, in complete agreement with (50), is that there must 
be rotational tori in this case. The median of the distribution of t c for 0.0512 < e < 0.10 can be fit 
to the power law 

t c = 0.0132(e - £oo) 



-3.807 



with = 0.05047 ± 0.0002 as shown in Fig. 19 The value is close to (50), and the difference 



is probably due to the weight given in the fit to points for larger values of e, though it may also be 
due to the sharp dependence of e cr on 8. 

Indeed, we compare the crossing time experiments for various values of 5 with the critical £ 
values of the tori in Fig. 19 'b). The points in the figure represent the computed (£ cr ,<5 cr ) for the 



2048 tori of Fig. 17 in the rrl and llr triangles. The upper envelope of these points is an estimate 
of the most robust torus for a fixed 5. Since these points are not from the zoomed- in triangles, like 
those of Fig. 18 they underestimate the maximal e cr . The crosses in Fig. [l~9^b) show the computed 
the positions of the pole, £oo, in the crossing time. These are obtained from power law fits to 
crossing times experiments with a maximal iteration time of 10 10 . Note that the maximal value of 
£ is a highly sensitive function of 5. 



The crossing time experiments confirm that the most robust torus of (16) for the standard 



parameters (18) is destroyed near (50). 




(b)- 




B-E K 10"" 



Figure 19: (a) Crossing time as a function of e for orbits of ( 16 1 with S = —0.33. Points correspond to the median 



of the crossing times for each e, and error bars to the median absolute deviation, (b) Critical e values as a function of 
8. The dots (black) are (e c 
Eoo, of t c for twenty values of 8. 



--cr, S cr ) pairs for the tori of Fig. |17| and the crosses (red) are estimates of the pole position, 
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9 Conjectures and Questions 



Our investigation of the reversible, two-angle, one-action map (16) has given support for several 
conjectures: 

1. For each rational rotation vector m/n in the range of £l(z, 5) and for each e, there is a 5 such 



that (16) has type a (m, n)-periodic orbit. This conjecture is the analogue of the Poincare- 



Birkhoff theorem, and seems to be open for the volume-preserving case. 
2. The residue (30) of each (m,n) orbit of (16) with force (15) grows as e or e 2 for eCl and 



as e n 1 or e n for e 3> 1. 

3. If the residue of a sequence of symmetric (m, n)-orbits limits to zero when m/n —> u for a 
Diophantine rotation vector uj, then there is a rotational invariant torus with that rotation 
vector. In other words, Greene's residue criterion holds. This conjecture is analogous to 
Tompaidis's theorem for the symplectic case |Tom 96a] . 

4. Conversely, if the residues of such a sequence are unbounded then the corresponding torus 
does not exist. This conjecture should be valid if the orbits converge to a remnant torus with 
positive Lyapunov exponent as shown by [FdlL92]. 

There are also many questions: 



Is there an analogue of the dominant symmetry line for reversible maps like (16)? More 
generally, is there a pattern to the signs of the residues and their oscillations for the various 
symmetric families of orbits? 

Does the critical spiral torus have a particular threshold residue analogous to Rth ~ 0.25 for 
the noble tori of area-preserving maps? 

Is there some significance to the oscillations in sign of the residues of the approximating 



periodic orbits, particularly near the cusps in the critical set shown in Fig. 16? 



Do the periodic orbits that are best approximants to a Diophantine rotation vector have any 
special significance in approximating the corresponding invariant torus? 

Are there remnant tori analogous to the cantori of twist maps? The remnant torus shown in 



Fig. 15 appears to be similar to a Sierpinski carpet. Is this true? It is known that symplectic 
maps of the form have invariant Cantor sets when £<1 for every rotation vector |MM92] . 
but no similar results are known for volume-preserving maps. 

Is there a cubic field that, like the golden mean for twist maps, gives rise to locally most 
robust tori in the volume-preserving case? Several other fields, in addition to the spiral field 
that we have studied, have been proposed |Loc921 lTom96b| ICFL04] . 
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Appendices 



A Reversors 

Suppose that a map / : M -)• M, with M = T d x M fc is reversible with involutions Si and S2, and 
let F be a lift of / to R d x M fc . 

Note that if S is a reversor for /, then so is S o /* for any ( 6 Z; indeed this follows because 
/* is a symmetry of /. For the lift F, the translations (20) are also symmetries; consequently, 
the translation of a reversor S is also a reversor, S o T m . Note that Si o T m = T_ m o Sj. These 
transformations generate a group 

(Si, S2, F, T m ) 



that is contained in the reversing symmetry group of F whenever g is odd. For (16) with (12) and 
(15), the reversors were given in (25), This appears, as far as we know, to be the complete group. 
The fixed set of S is 

Fix(S) = {(x, z) : S(x, z) = (x, z)} (51) 
As is well-known, every symmetric periodic orbit has points on the fixed sets. 

Lemma 4. If f has reversor S, then every symmetric orbit must have a point on either Fix(S) or 
Fix(/oS). 

Proof. Indeed, if j = 2k is even then 

(x k , z k ) = f k (x 2k , Z2k) = F k S(x , Zq) = S(x k , z k ), 
so that (x k , z k ) £ Fix(S). If, on the other hand, j = 2k — 1 is odd, then 

(x k , Zk) = f~ k+1 ° S(x , Zq) = f o S(x k , Z k ), 

so that (x k , z k ) G Fix(/ o S). □ 

Even better, periodic orbits have points on two distinct symmetry curves separated by half of 
the period. 

Lemma 5. If f : T d x M fe has a reversor S and T is a symmetric, (m,n) -periodic orbit of the lift 
F to the universal cover M. d x ~R k , then T has points on two fixed sets of reversors in the group 
generated by (S, F, T m ) . 

Proof. Since F is symmetric, we can suppose that there is a point (xq, zq) £ V (1 Fix(S). When the 
period n = 21 is even then, since (x_£, z_i) = T_ m (x£, zg), 

(xe, zi) = F e (S(x , z )) = S(x^i, z_i) = So T- m (x e , z e ); 

thus (x£,Zi) € Fix(SoT_ m ). Similarly when n = 2t — 1 is odd, then since (a;_^ + i, Z-i+i) = 
T- m (x£,Z£), 

(x e , ze) = F l (S(x , zq)) = F o S(x-e +1 , z-e+i) = F o S o T^ m (x e , z e ); 
thus (xe, zi) G Fix(F o S o T_ m ). □ 
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If we denote the symmetry by Si and let S2 = f ° S\, then, the lemmas imply that symmetric 
orbits can be found by looking for orbits that start on Fix(S'i) or Fix(S*2) = Fix(/ o Si), and "half 
a period later end on one of the sets Fix(5i o T_ m ) or Fix(S < 2 o T_ m ). The results are summarized 
in Table [Q 
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